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10 present application. The disclosure of each of these related applications is hereby 

7rporated by reference. 
Technical Field 

The present invention generally relates to the field of diagnostic spectroscopy, and 
more specifically, to an improved illumination device for obtaining spectral information 
15 and determining physical properties of a tissue sample. The present invention provides 
for an illumination device including an element or elements that spatially and angularly 
homogenize radiation for illuminating a sample under analysis while maintaining a high 
signal-to-noise ratio. 

Background of the Invention 
20 Spectral data arising from spectroscopic analysis provides practitioners with a 

wealth of detailed information about the identity, structure, and concentration of samples 
or constituents of samples. Spectral data derives from the detected and recorded energy 
change of a molecule through the emission or absorption of a photon. 



In particular, practitioners focus upon a molecule's vibration. Atoms within a 
molecular species vibrate back and forth about an average distance. Absorption of light 
by an atom at an appropriate energy causes the atoms to become excited, elevating the 
atom to a higher vibration level. The excitation of the atoms to an excited state occurs 
only at certain discrete energy levels, which are characteristic for that particular 
molecule. Infrared absorption spectroscopy is particularly useful for performing this type 
of analysis. In absorption spectroscopy, the net absorption of incident radiation at 
various wavelengths is measured. 

Radiation passing through a sample is attenuated depending upon the pathlength 
traveled by the radiation and the strength of absorptions at various individual 
wavelengths for constituents within that particular sample. Recording and mapping the 
relative strength of the absorption versus wavelength results in a unique absorption 
"fingerprint" for that particular sample. 

Cataloging infrared absorption permits practitioners to readily identify unknown 
samples by cross-referencing their samples of interest with a cataloged database. 
Matches between the spectrum of a compound of known structure and that of an 
unknown structure unambiguously identify the latter. This is an illustration of the 
qualitative aspects of spectrometry. Similarly, spectrometry also aids practitioners in 
quantitative analysis of known compounds. In illustration, calibration of a 
spectrophotometer on a known analyte of known concentration permits the accurate 
measurement of the same analyte of unknown concentration. 

Most absorption spectroscopic instruments consist of five components: 1) a 
source of radiant energy, 2) a wavelength selector that permits the isolation of a restricted 



wavelength region (i.e., a monochromator), 3) one or more sample containers, 4) a 
radiation detector which converts radiant energy to a measurable signal (usually 
electrical), and 5) a signal processor and readout. 

One application area for multivariate quantitative spectroscopy is the measurement 
5 of tissue attributes or analytes noninvasively. A specific application is the measurement of 
glucose noninvasively for subjects with diabetes. This application is difficult due to the 
complexity of the tissue, a turbid media, and the small size of the glucose signal. For the 
measurement of analytes with small concentrations in turbid media, care must be taken to 
minimize spectroscopic variances that overlap with the absorbance spectrum of the analyte 

10 of interest. Spectroscopic interferences are those spectroscopic variances unrelated to the 
analyte of interest but present during calibration development or during the measurement. 
Spectral interferences that overlap with or appear similar to the analyte of interest are 
especially bothersome. Spectroscopic interferences or variances can be classified into four 
general groups: 1) measurement variance; 2) physiological variance; 3) instrument 

15 variance; and 4) sampling variance. Measurement variance is the variance inherent in the 
data acquisition process. Johnson noise, shot noise, electronic noise and quantization errors 
are all types of measurement variance or measurement noise. Physiological variances are 
typically defined by the chemical or structural complexity of the sample. In the preferred 
application of noninvasive glucose measurement, significant physiological variance is 

20 present due to the complex nature of skin. Instrument variance is any variance due to 
changes in the performance of the optical measurement system. Changes in the 
performance of the illumination system would be a type of instrument variation. Sampling 
variance is due to errors associated with optical interfacing to the sample of interest. An 



-3- 



objective in designing optical measurement instrumentation is to maximize the net analyte 
signal. The net analyte signal is that portion of the pure component that is orthogonal to the 
other sources of spectroscopic variation. The pure component spectrum is the absorbance 
spectrum of the analyte of interest in the absence of other absorbing species. See, for 
example, "Net Analyte Signal Calculation in Multivariate Calibration," by Avraham Lorber, 
Analytical Chemistry Analytical Chemistry , Vol. 69, No. 8, April 15, 1997. 

Practitioners skilled in the art have designed measurement systems that minimize 
measurement variances associated with the recoding of optical information. In particular, 
Johnson noise and electronic noise may be minimized through effective instrument 
design. Additionally, the operation of modern A/D converters using 16-bit or higher 
digitization can substantially reduce the effect of digitization error. Thus, through 
instrument design, shot noise can be isolated as the predominant noise source in most 
modern spectrophotometer systems working in the near-infrared region of the 
electromagnetic spectra. 

Physiological variances are due to compositional or structural differences that 
vary in the biological sample but are unrelated to the analyte or attribute of interest. 
Significant complexity is present when the sample of interest is tissue. Biological tissue 
is commonly characterized as a turbid medium. Turbid media generally fail to permit 
any single ray of light from following an undisturbed pathway through the medium. In 
effect, turbid media are non-deterministic. That is, a light ray that enters a turbid medium 
may undergo several scattering events before finally exiting the medium. When many 
light rays are directed into a turbid medium, each of the exiting rays collected at any 
given point will have traveled a different distance through the medium. As a result, a 



spectrum from a turbid medium source is determined not only by type and concentration 
of the molecular species in the medium, but also by the scattering characteristics of the 
media that influence the path distribution of the light traveling through the medium. 

Instrumentation variances are spectroscopic variations introduced by the optical 
5 instrumentation. Instrumentation variances can result in statistically relevant errors. A 
recognized instrumentation variant is the illumination system. The illumination system is 
the entire optical system from the source of radiation to the sample or the spectrometer. 
Common components of an illumination system include the emitter or radiation source, 
reflective optics, refractive optics, focusing lenses, collimating lenses, filters, relay optics 

10 or optical fibers. A radiation source for absorption measurements must generate a beam 
with sufficient power in the wavelength region of interest to permit ready detection and 
measurement. In addition, the radiation source must also provide a reproducible output. 

The most common source of near-infrared radiation is the tungsten-filament lamp. 
The energy distribution of this source approximates that of a blackbody, and is thus 

15 temperature dependent. In most absorption spectrophotometers, the operating filament 
temperature is about 2900 K; therefore, the bulk of the energy is emitted in the infrared 
region. A tungsten-filament lamp is useful for the wavelength regions between 320 and 
2500 nm, the lower limit imposed by absorption by the glass envelope that houses the 
filament. 

20 A change in the lamp, adjustment to the lamp or changes within a lamp may affect 

the resulting spectral data produced by the illumination system. The following are four 
examples of common variants associated with the illumination system. The replacement 
of the lamp can result in significant spectroscopic variance when using existing 
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illumination systems. The spectroscopic variance can be due to manufacturing 
inconsistencies between lamp filaments. Known manufacturing inconsistencies include 
differences in filament shape, differences in filament location, and differences in filament 
material. A second source of spectroscopic variance can be caused by rotation or tilting 
of the same lamp in the lamp housing. A third known source of variance is due to 
differences in the glass envelope surrounding the filament. Specifically the glass 
envelope "nipple" can create shadowing and cause inhomogenous illumination of the 
sample or sampling apparatus. Finally, the lamp may change over time due to vibration 
or sagging of the filament. Such changes can cause intensity and temperature variations 
along the filament length. With current illumination systems, the above changes can 
cause spectroscopic variances. In maximizing overall system performance, it is desirable 
to minimize spectroscopic variances unrelated to the analyte of interest. Illumination 
system variances include all variances due to different lamps, due to lamp aging, due to 
placement of the lamp in the instrument or any other variance that results due to a change 
in lamp performance or how the radiation source interacts with the remainder of the 
optical system. With current illumination systems, radiation emitter variances can cause 
spectroscopic variances. Illumination system variances unrelated to either the sample of 
interest or the analyte being measured can result in prediction errors and necessitate the 
need for re-calibration. Re-calibration is generally undesirable due to increased expense 
and down time on the instrument. Thus, in maximizing overall system performance, it is 
desirable to minimize illumination system variances. 

To achieve increased accuracy in the measurement of analyte concentration, a 
practitioner of the art must, among other things, strive either to eliminate interferents or to 



construct a chemometric model that is sensitive to the differences between interferents and 
the desired analyte. Fortunately, spectral changes due to interferents are seldom identical to 
spectral changes due to changes in analyte concentration. Thus, the ability of a chemometric 
model to distinguish between an interferent and a particular analyte is typically improved by 
increasing the size of the calibration set. 

hi the presence of significant spectroscopic variance and when the degree of overlap 
(spectral similarity) between the analyte and the interferent is high, the number of model 
"factors" required to adequately distinguish between the interferent and the analyte will be 
large (the model complexity will be high). Unfortunately, in the presence of measurement 
noise, there are practical limitations associated with the number of model factors that can be 
used effectively. The ability of the practitioner to improve the sensitivity of the model to 
differences between the interferent and the analyte by increasing the size of the calibration 
model will be limited by the presence of noise in the measurement which limits the ability to 
distinguish between the spectra of the interferent and the analyte. In systems where there is 
a high degree of overlap between interferents and the particular analyte of interest, the 
practitioner must strive to reduce the amount of spectroscopic interference, specifically 
instrument variance to the greatest degree possible. 

In spectrophotometer instruments where shot noise is the predominant source of 
measurement noise in the instrument, the signal-to-noise ratio (SNR) for the instrument is 
directly proportional to the square root of the flux (0 ) on the photodetector. Thus, for these 
instruments, the SNR can be improved by maximizing the amount of light incident on the 
detector. For measurements on biological tissue, however, the practitioner cannot increase 
the flux on the detector without limit. Increasing the flux on the detector generally 
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necessitates increasing the incidance on the tissue. The increased incidence on the tissue 
may result in thermal damage to the tissue. Therefore, there are practical limits on how 
much light can be incident on the tissue. 

Fourier-transform infrared (FTIR) spectrophotometers are a class of spectrometer 
5 that can be operated where shot noise is the predominant measurement noise. FTIR 
spectrophotometers offer the advantages of unusually high sensitivity, resolution, and speed 
of data acquisition. Generally, data from an entire spectrum can be obtained in one second 
or less. The heart of a Fourier transform spectrophotometer is an interferometer, which is a 
device for analyzing the frequencies present in a composite signal and the relative strength 
10 of the signal at such frequencies. 

Vast improvement in spectroscopic optical analysis, particularly quantitative 
determination of analytes in biological tissue, can be achieved if the above-identified 
problems are minimized or eliminated. In particular, if illumination system variances 
could be eliminated or substantially reduced as an interferent the complexity (number of 
15 model factors) of the model may be reduced and the net analyte signal increased. This 
would result in reduced instrument variance and increase the predictive accuracy of a 
model, particularly in situations where the interferent has substantial overlap with the 
analyte of interest in the tissue. It is further important that a system which eliminates or 
reduces instrument interferents also minimizes measurement noise. The present 
20 invention is directed to apparatus and methods which eliminate or reduce changes in the 
light source or illumination system as an interferent while maintaining high signal-to- 
noise ratio. 
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Summary of the Invention 
The present invention provides an apparatus and method for spatially and 
angularly homogenizing radiation emitted by a radiant source for use in illuminating a 
sample or sampling apparatus for spectroscopic optical analysis. More specifically, the 
5 present invention provides components for spatially and angularly homogenizing 
radiation emitted by a radiant source which is to be used to illuminate a sample under 
analysis that overcomes expected variances particular to that radiant source and its 
placement and use within an optical system. Further, the present invention provides an 
apparatus and method for improving signal-to-noise ratio by filtering selected 
^ 10 wavelengths prior to contact with the sample. 

gj The reduction of the illumination source as an instrument variance or interferent 

fn has been found to greatly improve the ability to build an optical system and model which 

m 

can accurately predict small analyte concentrations in turbid media such as tissue. The 
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present invention provides this illumination stability not by use of a new radiation source, 
15 but by modifying the output beam prior to sample or sampler illumination to minimize 
spectroscopic variances due to the illumination system. 

The present invention relates to an apparatus for minimizing spectroscopic 
variances due to radiation emitter differences by use of an illumination system that 
utilizes both angular and spatial homogenization. Angular homogenization is any process 
20 that takes an arbitrary angular distribution, or intensity (W/sr), of emitted radiation, and 
creates a more uniform angular distribution. Spatial homogenization is the process of 
creating a more uniform distribution of irradiance (W/m 2 ) across an output or exit face. 
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A standard lamp produces a non-uniform irradiance distribution due to the 
physical structure of the lamp filament. Thus, radiation emitter differences (e.g., a 
different lamp) will result in different non-uniform irradiance distributions. These 
differences in irradiance distribution between the lamps can translate into spectroscopic 
5 differences. Thus, an objective of the invention is to take different irradiance 
distributions due to emitter differences and create similar or ideally the same irradiance 
distribution. A preferred method of creating similar irradiance distributions is to create a 
uniform irradiance distribution. 

Differences in the radiation emitter will also result in differences in angular 

10 distribution. As above, an objective of this invention is to create an illumination system 
where radiation emitter differences do not affect the angular distribution observed by the 
sample or at the input to the spectrometer. One mechanism is to create a uniform angular 
distribution. An ideal angular homogenizer would uniformly distribute the light over a 
sphere (4 pi sr) regardless of the angular distribution from the emitter. An ideal reflective 

15 angular homogenizer would uniformly distribute light over a hemisphere (2 pi sr). Due 
to the fact that other optical components in the system must collect light within a defined 
numerical aperture, ideal diffiisers are typically very inefficient. Thus, the instrument 
designer must weigh the benefits of angular homogenization with loss in optical 
efficiency. Regardless of the specific embodiment, angular homogenization is a critical 

20 component in the realization of an illumination system that has reduced sensitivity to 
emitter differences. 

The present invention provides a system for producing spatially and angularly 
homogenized light from an irregular emitter and using the homogenized light for spectral 
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analysis. The resulting homogenized radiation illuminates the sample or sampler in a 
consistent and reproducible form, thus allowing for accurate and dependable molecular 
absorbance measurements. 

An additional benefit of the current invention is spectral or "color" 
5 homogenization. A typical quartz tungsten halogen bulb emits blackbody radiation in a 
fashion that can be described by Planck's radiation law. This emitted radiation is a 
function of the filament's color temperature, a colorimetric concept related to the 
apparent visual color of a source. The color temperature of the filament is not constant 
across the length of the filament and therefore the spectral radiation will not be constant. 

10 Thus, color temperature variations across the filament will result in spectral differences 
across the filament length. These spectral differences due to color temperature variations 
or other filament differences can be different between emitters and can change over time. 
These differences in color temperature or spectral distribution between lamps can 
translate into spectroscopic differences. Thus, an additional objective of the invention is 

15 to take different spectral distribution due to emitter differences and create similar or 
ideally the same spectral distribution. A preferred method of creating similar spectral 
distribution is to create an output that has uniform spectral content. 

The usefulness of the present invention is best illustrated by the familiar 
occurrence of routine maintenance to a spectrometer. It is common for radiant light 

20 sources to burn out. Although application dependent, the replacement of the light source 
has resulted in prediction errors and has necessitated recalibration of the spectrometer to 
the new light source. In systems intended for commercial use by unskilled operators, 
recalibration is not desired. With the present invention, however, differences in the light 
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source are irrelevant and proper performance of the optical measurement system is 
maintained. Regardless of the spatial and angular characteristics of the radiation entering 
the light homogenizes the use of the present invention will result in radiation incident on 
the sample which always remains generally spatially and angularly homogenized. Thus, 
a light source change will not detract from the accuracy and dependability of molecular 
absorbance measurements using the present invention. 

The present invention further specifies a system for providing illumination to 
biological tissue samples. More specifically, the system is particularly suited for 
spectroscopic illumination of biological tissues for determining and quantifying the 
concentration of specific analytes within or other characteristics of the tissue. The present 
invention enables a practitioner to construct and operate an illumination device that permits 
measurements with a high signal-to-noise ratio (SNR) while minimizing thermal damage to 
biological tissue. With a high SNR, chemometric models may be developed for 
differentiating between a particular analyte and interferents similar to that analyte. The 
present invention allows for infrared spectroscopic analysis of turbid medias by satisfying 
the following conditions: 

(1) The radiation emitted by the present invention contains wavelengths useful for 
measuring the analyte of interest. The radiation may be continuous versus wavelength, in 
locally continuous bands, or selected to particular wavelengths. The result is radiation that 
encompasses the wavelength regions that contain the near-infrared spectral "fingerprint" for 
the analyte of interest. For FTIR spectroscopy of glucose, this wavelength region spans 
approximately from 1.4 to 2.5 (am. 
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(2) The radiation emitted by the present invention is of sufficiently high spectral 
radiance (L x ) to provide a high signal-to-noise ratio in the spectral fingerprint region for the 

analyte of interest. In FTIR spectroscopy of glucose, for example, the radiation from a 
halogen lamp concentrated with one or more optical elements, such as lenses and or mirrors, 
will provide a spectral radiance that satisfies this condition. 

(3) The spectral incidance (E x ) on the tissue is low enough to avoid burning 

biological tissue, or otherwise (discomforting) subjects. This may include the use of spectral 
filters to reduce or eliminate radiation that is outside the fingerprint region, and to 
selectively suppress certain wavelength bands that lie within the fingerprint region. Spectral 
filtering may be achieved by absorptive, diffractive, or interference means (absorption 
filters, gratings or prisms, or multilayer dielectric filters). 

(4) The spectral radiance (L x ) is generally invariant when subjected to changes in 

the spectral exitance (M \) of the emitter. Reasonably expected changes in the spectral 

exitance are those due to rotation and/or small translation of the emitter, or replacement of 
the emitter with another emitter of the same general construction. 

By satisfying the above conditions, the present invention eliminates the need for 
recalibration due to illumination variability (bulb changes, source rotation or movement) 
or development of a chemometric model that compensates for such changes. Simple 
maintenance such as replacing the light source would no longer necessitate recalibration 
or the development of chemometric models sensitive to light source changes. 
Furthermore, rotations and translations of the light source caused by jolts, bumps, and 
other similar vibrations would have minimal effects on the accuracy of a test. 



-13- 



Figure 1 is a detailed perspective view of an infrared radiation source lamp known 
in the art; 

Figure 2 is a diagrammed view of a system for measuring the concentration of an 
analyte within biological tissue; 

Figure 3a is an incidance plot using a ray trace program simulating the spatial 
distribution of emitted radiation from an infrared spectrophotometer known in the art; 

Figure 3b is an incidance plot showing the changes in spatial distribution of 
emitted radiation after a 90-degree rotation of the filament used in producing the 
incidance plot of Figure 3a; 

Figure 3c is an incidance plot showing the changes in spatial distribution of 
emitted radiation after a one-millimeter vertical translation of the filament used in 
producing the incidance plot of Figure 3 a; 

Figure 4a is an intensity plot using a ray trace program simulating the angular 
distribution of emitted radiation from an infrared spectrophotometer known in the art; 

Figure 4b is an intensity plot showing the changes in angular distribution of 
emitted radiation after a 90-degree rotation of the filament used in producing the intensity 
plot of Figure 4a; 

Figure 4c is an intensity plot showing the changes in angular distribution of 
emitted radiation after a one-millimeter vertical translation of the filament used in 
producing the intensity plot of Figure 4a; 

Figure 5 is a diagrammed view of a system for constructing a chemometric model 
for measuring glucose concentration in the forearm's of various subjects; 
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Figure 6 is a box and whisker plot of prediction error versus day across five lamp 
changes using the system illustrated in Figure 5; 

Figure 7 is a box and whisker plot of in-vivo prediction errors versus orientation 
for a lamp within a system illustrated in Figure 5; 

Figure 8 is a diagrammed view of a system used for cross-validation analysis for 
baseline system performance using a tissue phantom for the sample source; 

Figure 9a is a box and whisker plot of cross-validated prediction errors for the 
system illustrated in Figure 8, in the absence of a lamp change; 

Figure 9b is a box and whisker plot of cross-validated prediction errors for the 
system illustrated in Figure 8, with the inclusion of lamp changes; 

Figure 10 is a diagrammed view of a system of the present invention using a 
means for spatially and angularly homogenizing emitted radiation; 

Figure 11a and Figure 1 lb illustrate a detailed perspective and plan view of a light 
pipe of the present invention; 

Figure 12 is a plan view of a ray trace showing radiation focused by an elliptical 
reflector into and through a light pipe of the present invention; 

Figure 13a is an incidance plot using a ray trace program simulating the spatial 
distribution of emitted radiation using a light pipe of the present invention; 

Figure 13b is an incidance plot showing the changes in spatial distribution of 
emitted radiation after a 90-degree rotation of the filament used in producing the 
incidance plot of Figure 13 a; 



-15- 



Figure 13c is an incidance plot showing the changes in spatial distribution of 
emitted radiation after a one-millimeter vertical translation of the filament used in 
producing the incidance plot of Figure 13a; 

Figure 14a is an intensity plot using a ray trace program simulating the angular 
distribution of emitted radiation using a light pipe of the present invention; 

Figure 14b is an intensity plot showing the changes in angular distribution of 
emitted radiation after a 90-degree rotation of the filament used in producing the intensity 
plot of Figure 14a; 

Figure 14c is an intensity plot showing the changes in angular distribution of 
emitted radiation after a one-millimeter vertical translation of the filament used in 
producing the intensity plot of Figure 14a; 

Figure 15 is a schematic plan view of an alternative source and light pipe system 
of the present invention; 

Figure 16 is an incidance plot depicting homogenization of the light at the distal 
end of the light pipe of Figure 15; 

Figure 17 is an intensity plot showing the homogenization of light emitted from 
the light pipe of Figure 15; 

Figure 18 is a schematic plan view of an alternative illumination source 
incorporating parabolic reflectors and a light pipe; 

Figure 19 is an incidance plot depicting spatial homogenization of the light; 

Figure 20 is a plot of intensity showing the homogenization of light by the source 
in Figure 18; 
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Figure 21 is a schematic perspective view of an alternative illumination source 
incorporating faceted reflectors; 

Figure 22 depicts spatial distribution of the light showing spatial homogenization 
achieved through the system of Figure 21; 

Figure 23 is a plot of angular distribution produced by the device of Figure 21 ; 

Figure 24 is a diagrammed view of a system of the present invention for 
measuring glucose in scattering media having a tissue phantom as the sample source; 

Figure 25a is a box and whisker plot of a standard system with no bulb changes; 

Figure 25b is a box and whisker plot of a standard system across four bulb 
changes; 

Figure 25c is a box and whisker plot of a system using an s-bend light pipe across 
four bulb changes; 

Figure 25d is a box and whisker plot of a system using a ground glass diffuser 
plus and s-bend light pipe across four bulb changes; 

Figure 26 is a diagrammed view of a system incorporating filters prior to the light 
pipe which eliminate unwanted wavelengths of radiation from the illumination source; 
and 

Figure 27 graphically depicts the transmittance of selected wavelengths in a 
preferred fingerprint region. 

Detailed Description of the Preferred Embodiments 

The following detailed description should be read with reference to the drawings, 
in which like elements in different drawings are numbered identically. The drawings, 
which are not necessarily to scale, depict selected embodiments and are not intended to 
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limit the scope of the invention. Examples of constructions, materials, dimensions, and 
manufacturing processes are provided for selected elements. Those skilled in the art will 
recognize that many of the examples provided have suitable alternatives that may be 
utilized. 

5 Figure 1 shows a plan view of a near infrared radiation source lamp 14 known in 

the art. The appearance of a radiant source lamp 14 closely resembles that of a traditional 
residential light bulb. Traditional spectrophotometer lamps consist of a filament 16 
housed within a transparent envelope 18, or the like. The transparent envelope 18 is 
either comprised of a silicate glass, fused silica or quartz material. The material used for 

10 the glass envelope 18 is dependent upon the wavelength regions being surveyed on the 
electromagnetic spectrum. 

The envelope 18 traditionally is cylindrical or oval in shape. The lamp 14 of 
Figure 1 specifically is of a closed-end cylindrical variety. The closed-end portion of the 
cylinder has a nipple 20 positioned near the center of the cylinder's closed-end base. The 

15 nipple 20 formation is a result of manufacturing and functionally has no beneficial 
purpose. On the other hand, the nipple 20, as will be discussed in detail later, affects the 
emission of radiant energy. 

Filament 16, and subsequently lamp choice, is wavelength dependent. Operating 
in the infrared and near infrared regions of the electromagnetic spectrum requires a 

20 radiation source filament 16 applicable to those spectral regions. Several continuous 
radiation sources including tungsten-halogen lamps, tungsten lamps, nerst glowers, 
nichrome wires and globars are suitable for infrared molecular absorption spectroscopy. 
The desired filament is manufactured so as to place the filament 16 within the open end 
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of the glass envelope 18 and securely fastened thereto. Wires or leads 22 emerge from 
the filament 16 and out of the glass envelope 18 attaching the filament 16 to an energy 
source (not shown). Because the energy output of a filament 16 generally varies 
approximately with the operating voltage, close voltage control is essential. For this 

5 reason, most lamps 14 are attached through the wires or leads 22 to a constant-voltage 
transformer or electronic voltage regulator. 

Tungsten-halogen lamps are a preferred radiation source for performing infrared 
molecular absorption spectroscopy. Tungsten-halogen lamps are particularly useful in 
modern spectroscopic instruments because of their extended wavelength range, greater 

10 intensity, and longer life. Specifically, tungsten-halogen lamps emit wavelengths within 
the region of 240-2500 nm, covering portions of the ultraviolet, visible, and near-infrared 
spectral regions of the electromagnetic spectrum. Thus, a tungsten-halogen lamp is a 
preferred source of radiant energy for the present invention. Any suitable lamp source 
operating within the infrared region of the electromagnetic spectrum, however, may be 

15 used in conjunction with the present invention. 

The basic illumination source depicted in Figure 1 further includes an elliptical 
reflector 12 which focuses emitted light from bulb 14 to a reflector focus 26. 
Representative rays 24 are depicted to show the function of the reflector 12. The 
relationship between the radiant source emitter 14 and the elliptical reflector 12 was used 

20 in the subsequently disclosed experiments. 

Referring now to Figure 2, a diagrammed view of a system 10 for measuring the 
concentration of an analyte within biological tissue is depicted. The system 10 depicted 
is simplified by illustrating certain specific elements within a far more elaborate 
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spectroscopic system. The elements depicted in Figure 2, however, are common to 
spectroscopic systems, and therefore, require some identification. 

An elliptical reflector 12 known in the art is shown. At the center of the elliptical 
reflector 12 is radiation source lamp 14. The radiation source lamp 14 is depicted as 
having a filament 16, a glass envelope 18 with nipple 20 housing the filament 16, and a 
pair of leads 22 extending from the end of the lamp. Surrounding a portion of the lamp 
14 is the body of the reflector 12. The elliptical reflector 12 functions to concentrate 
emitted radiation rays 24 (shown as a ray trace) from the radiation source lamp 14 onto 
the reflector's focal point 26, In order to maximize reflectance, the elliptical reflector 12 
is generally made from a highly polished metal. Although Figure 2 specifically illustrates 
an elliptical reflector, other shapes suitable for focusing radiant energy are also within the 
scope of the invention. 

Figure 2 depicts two fiber optic bundles, an illumination fiber bundle 30 and a 
collection fiber bundle 32. Fiber optic bundles 30 and 32 are extremely versatile because 
they are capable of channeling harnessed radiation between elements without noticeable 
reduction in the intensity of that harnessed radiation. 

At the reflector's focal point 26 is an opening to the illumination fiber bundle 30. 
The illumination fiber bundle 30 collects the radiation emitted 24 by the lamp 14 and 
channels the radiation through the bundle system. At the output end of the illumination 
fiber bundle 30 is another opening that then directs the harnessed radiation onto a sample 
40, such as human tissue on a person's forearm. The second fiber optic bundle, the 
collection fiber bundle 32, is positioned proximate to the sample 40 to again collect 
radiation, however, here the radiation is diffusely reflected from the sample 40. 
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Diffusely reflected radiation is that radiation reflected from within the sample 40. 
Diffusely reflected radiation does not generally follow a uniform pattern. Ray tracing of 
the diffusely reflected radiation within the sample 40 as shown in Figure 2 illustrates 
possible pathways of radiation entering, and subsequently reflecting out of, the sample 
5 40. 

The sample 40 is the subject of the spectroscopic analysis. Qualitative and 
quantitative measurements of specific analytes within the sample 40 may be conducted 
by measuring absorption of introduced radiation from the emitter 14. In a preferred 
embodiment, the sample 40 is biological tissue. More specifically, the sample 40 is an 

10 appendage, or a portion thereof, of the human anatomy. For example, the sample 40 may 
be a forearm or a finger, to list a few. 

The collection fiber bundle 32 then channels the diffusely reflected radiation from 
the sample 40 to a spectrophotometer 44 where the radiation is detected by converting the 
recaptured radiant energy into a measurable signal. In a preferred embodiment of the 

15 invention, a Fourier transfer infrared (FTIR) spectrophotometer is utilized to analyze the 
diffusely reflected radiant energy emitted by the sample 40. The usually electrical signal 
is then transferred to a signal processor. Processing of the signal is generally 
accomplished using a computer or other data processing means 46 designed for such 
processing. The outcome of the processing is then transcribed to a readout, allowing 

20 practitioners to study the results of the analysis. 

As described in detail above, in spectrophotometer instruments where shot noise 
predominates the system, as is in the system depicted in Figure 2, the signal-to-noise ratio 
(SNR) for the system is directly proportional to the square root of the flux (0) on the 
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photodetector. The SNR for the system, however, can be improved by maximizing the 
amount of radiation incident on the detector. Increasing the flux on the detector generally 
necessitates increasing the incidance, and thus, may cause thermal damage on the 
sampled biological tissue 40. To illustrate this tissue-heating problem, experiments were 
conducted utilizing the system illustrated in Figure 2. For the experiment, the sample 40 
used was the forearm of a living human subject and the analyte to be measured was 
glucose. 

The radiation source lamp 14 was connected to a variable current source that 
permitted the lamp 14 to increase output up to a maximum of 40 watts. The output of the 
lamp would then be incrementally increased until the SNR was high enough to acquire 
accurate glucose measurements. As the lamp power was increased during the subsequent 
experimental trials, most of the subjects reported discomfort prior to reaching an 
acceptable SNR. The discomfort experienced by the subjects was due to a localized 
heating of their forearm by the illuminating radiation. 

To further analyze the above-described phenomenon, a ray trace program was 
utilized to compare and contrast various illumination systems for spatial and angular 
homogeneity. TracePro V2.1, a commercially available non-sequential ray trace 
program, was used to generate realistic models of the radiation distributions from various 
illumination system configurations. The output from such modeling is depicted in 
Figures 3a-c, 4a-c, 13a-c and 14a-c. In order to understand the output of the modeled 
illumination, Table 1 correlates the specific radiometric terms to their corresponding 
symbols, definitions, and units. 
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TABLE 1: Definition of Radiometric Quantities 



Name 


Symbol 


Definition 


Units 


Energy 


. Q 




Joules, J 


Flux 


0 


8Q 
dt 


Watts, W 


Exitance 


M 


80 
8A S 


W/m 2 


Incidance 


bi 


80 
8A r 


w/m 


Radiance 


L 


80 

d(A s cos0) dn 


W/m 2 /sr 



With respect to Table 1, dA r and dA s refer to differential elements of area on the 

receiver and source, respectively. Additionally, 6 refers to the angle between the line of 

sight from the observer to the source and the direction of the radiation 24. The associated 

spectral quantities are defined by differentiating the above general radiometric quantities 

with respect to wavelength, as depicted below: 

M dM dE dL 

M ; = , E 2 = — , and L } = — 

A dX k dX 1 dX 

Figures 3a-c are plots of the incidance of emitted radiation 24 from the elliptical 
reflector 12 in Figure 1. These plots have been generated using TracePro V2.1, a ray 
trace program simulating the spatial distribution of emitted radiation from the radiation 
source lamp 14. More specifically, the plots of incidance are representative of the spatial 
distribution of emitted radiation at the focus of the elliptical reflector 26 diagrammed in 
Figure 1. 
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Figure 3a shows a plot of incidance of emitted radiation 24 from a radiation 
source lamp 14. The resulting incidance plot is characterized by a substantial degree of 
spatial inhomogeneity. Spatial distribution of emitted radiation in particular areas of the 
plot is demonstrated to vary substantially throughout the incidance plot. In certain areas 
within the plot, the spatial distribution is greater than other areas within the same plot. 
The converse is also true. The spatial distribution of the emitted radiation is also 
illustrated to follow certain arc-like bands of greater or lesser incidance throughout the 
plot. 

Figure 3b shows a plot of incidance of the same radiation source lamp of Figure 
3a, but after a 90-degree rotation of the filament producing the incidance plot. 
Comparisons of the plots of Figures 3a and 3b show that areas of greater incidance in 
Figure 3a are now areas of lesser incidance in Figure 3b, and the inverse. Figure 3c 
further depicts this spatial distribution disparity by showing the changes in spatial 
distribution when the filament 16 of the same radiation source lamp 14 of Figure 3a 
undergoes a vertical translation of one millimeter. Again, the spatial distributions in 
Figure 3c after the one-millimeter translation provide areas of greater incidance where 
there were originally none in Figure 3a. These plots document that the spatial 
distribution of light at the focus of the standard light source is highly unstable with 
modest translations and/or rotations of the filament. 

Similar to Figures 3(a-c), Figures 4(a-c) depict plots of the intensity of emitted 
radiation from the elliptical reflector in Figure 1. These plots have also been generated 
using TracePro V2.1 to simulate the angular distribution of emitted radiation 24 from a 
radiation source emitter 14 known in the art. More specifically, the plots are 
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representative of the angular distribution of radiation at the focus of the elliptical reflector 
26 diagrammed in Figure 1, i.e., the direction of the light rays at the focus of the elliptical 
reflector. 

Figure 4a shows a plot of intensity of emitted radiation from a radiation source 
lamp 14. The resulting intensity plot from the standard radiation source is characterized 
by a substantial degree of angular inhomogeneity. Angular distributions in particular 
areas of the plot also vary dramatically within the same plot. For example, Figure 4a 
illustrates a "hole" in the center of the intensity plot. The lack of irradiation intensity in 
this particular area is a result of a shadowing effect by the envelope nipple 20 on the end 
of a radiation source lamp 14. 

Rotating the filament 16 of Figure 4a produces an intensity plot illustrated by 
Figure 4b. Because the filament 16 was rotated, the hole 60 in the center of the plot 
remains centered within the plot after the 90-degree rotation. Translation of the filament 
16 of Figure 4a by one millimeter, however, greatly diminishes the angular distribution 
within the spectroscopic system, as depicted in Figure 4c, Angular distributions are 
sporadic, and often completely shadowed by the modest translation of the radiation 
source lamp 14. 

The ray trace plots of Figures 3(a-c) and Figures 4(a-c) illustrate that the spatial and 
angular distribution of light at the focus 26 of a standard radiation source 14 is highly 
unstable with respect to modest translations and/or rotations of its filament 16. Areas of 
higher incidance and intensity may form "hot spots" during illumination. In an attempt to 
maximize the signal-to-noise ratio (SNR), the radiation source 14 could be increased to the 
thermal and/or comfort limits established by the patient. However, if there are "hot spots" 
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across the tissue, these areas may require a lower overall radiation output and corresponding 
result of lower SNR. Thus, uniform intensity illumination is desired when attempting to 
maximize the SNR for tissue measurements such as glucose. 

The above plots clearly illustrate angular and spatial variances associated with the 

5 illumination system. These variances translate into spectroscopic variances that adversely 
influence the achievement of high levels of accuracy in measuring analyte concentrations. 
Inhomogeneous spatial and angular distributions of emitted radiation 24 impede a 
practitioner from constructing chemometric models that are sensitive to the differences 
between interferents and the desired analyte. Modest and unaccounted for translations 

10 and/or rotations of the emitter 14, such as those that might result from loose mechanical 
tolerances or vibration, have been found to significantly alter these relied upon chemometric 
models. An additional experiment was conducted to illustrate the effect of interferent 
variations on a calibrated chemometric model. 

Figure 5 shows a diagrammed view of the system used for constructing a 

15 chemometric model for measuring glucose concentration in the forearm of various subjects. 
The components within this instrument system closely resemble those in Figure 2 and like 
elements are numbered the same. The additions utilized should not be construed as an 
exhaustive list for constructing an accurate chemometric model for glucose measurement. 
Identification of these additions is merely for illustrative purposes only, as one of skill in the 

20 art may readily identify numerous combinations of instrument components that could 
achieve a chemometric model for the desired analyte. 

The first of the additions shown in Figure 5 is a five (5) millimeter aperture 70 
positioned at the focal point of the elliptical reflector 26. This aperture 70 limits the amount 



-26- 



of emitted radiation 24 permitted to pass through the system 10 for analysis. Once the 
radiation clears the aperture 70, a silicon lens 72 redirects the radiation through a cyan filter 
76, which in turn, sends the radiation through a second silicon lens 74. Radiation 
transmitted through this series of lenses is then filtered to absorb radiation at wavelengths at 
or greater than 2.7 micron by passing through filter/diffiiser 78. In a preferred embodiment, 
a WG295 filter/diffiiser is utilized to absorb the wavelengths at or greater than 2.7 micron. 
The radiation is then illuminated upon a sample 40, collected and analyzed as described in 
relation to Figure 2. 

Using the above-described system, numerous calibration spectra spanning a 
wavelength range of approximately 1.25 |im to 2.5 \xm were used to construct a 
chemometric model for measuring glucose concentrations within forearms 5 of subjects. The 
calibration set spanned several different lamps, many human subjects, a wide range of 
glucose values, and a variety of operating temperatures, and relative humidities. 

During the "prediction" phase of the experiment, eleven human subjects were 
measured by the spectrometer system four times each day. Additionally, the radiation 
source lamp 14 for the system was changed every two days. As a note, the human subjects 
and lamps used in this prediction phase of the experiment were not the same as those used 
during the calibration phase. The results of this experiment are shown in Figure 6, where 
the errors are sorted by day. 

Figure 6 shows a "box and whisker" plot. In this type of plot, the median prediction 
error for each day is plotted as a horizontal line 82 in the middle of a box 80, which 
encompasses the middle half of the data, and 'Svhiskers" 84 are plotted at the 5th and 95th 
percentiles; a "dot" 86 represents the mean prediction error for the day; the horizontal 



-27- 



dashed line 88 shows where the data are centered when the prediction error bias is zero; and 
the numbers shown at the bottom of the graph indicate the number of predictions associated 
with that whisker and taken on each study day. 

Figure 6 specifically shows a box and whisker plot of prediction error versus day 
across five lamp changes, six lamps in total, over twelve days. During the first four days of 
the experiment, regarding lamps 1 and 2, the absolute prediction error bias was less than 
20 mg/dl. After the second lamp change, however, (on days 5 and 6 of the experiment) the 
absolute bias increased dramatically. Replacing the third lamp with a fourth (on day 7 of the 
experiment) reduced the bias to well under 20 mg/dl. 

These results suggest that the chemometric model was sufficiently "robust" as to 
permit accurate determination of the glucose levels for the subjects for most of the lamps, 
even though the lamps used during the prediction phase were not the same as those used 
during calibration. With regard to the third lamp, however, the chemometric model failed to 
produce accurate predictions. This failure suggests that the emission characteristics of this 
lamp were substantially out of the calibration range used to build this experimental 
chemometric model. 

To help isolate the emitter variation as the source of the prediction errors described 
above, another experiment was conducted using the same apparatus, and similar methods as 
described in the previous experiment. In this subsequent experiment, however, spectra were 
collected from three different subjects all on the same day, using the same lamp throughout 
the prediction period. The lamp was installed in the apparatus at some arbitrary azimuthal 
orientation, 6 0 , and spectra of the subject's forearms were taken at 0 0 , as well as at 6 0 +/- 
2 degrees. The resulting prediction errors are plotted in Figure 7 for the three lamp 
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orientation states. These results indicate that changes in the emitter characteristics, which 
are the result of small rotations of the lamp, can cause prediction errors that are almost as 
large as those caused by complete replacement of one lamp with another. 

A third experiment was then conducted to evaluate the effects of lamp changes on 
prediction error. The system utilized is depicted in Figure 8, with like elements numbered 
the same as in Figure 2 and Figure 5. In this experiment, the sample source of living tissue 
40 (a subject's forearm) was replaced with a "tissue phantom" 43. Tissue phantoms 43 
consist of a scattering solution made of microscopic polystyrene beads suspended in water at 
varying concentrations. In this experiment, the concentration range for the polystyrene 
beads was between 5000-8000 mg/dl. Tissue phantoms 43 within these ranges are 
representative substitutes for living tissue because their optical scattering and absorption 
properties are similar to those of biological tissue. Additionally, the use of tissue phantoms 
of known concentrations eliminates the confounding effects often observed from 
physiological changes in living tissue. Figure 8 diagrams the replacement of a subject's 
forearm 40 with a tissue phantom 41. Further, the cyan filter 76 is located after the output 
fiber optic 32. In all other respects, the apparatus diagrammed in Figure 8 is consistent with 
those discussed in detail with respect to Figure 2 and Figure 5. 

A set of ninety-eight different tissue phantoms composed of five different analytes at 
different concentrations was optically sampled. In order to assess the ability of the system 
in Figure 8 to predict glucose concentrations in the absence of lamp changes, a "cross- 
validation" analysis was performed. To accomplish this cross-validation analysis, a series of 
baseline measurements were performed wherein spectra of all ninety-eight solutions were 
taken using a single lamp with the apparatus depicted in Figure 8. This data was artificially 
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subdivided into four sets. Using three of these sets, a chemometric model was constructed 
to predict glucose values for the remaining set. The analysis procedure was again repeated, 
rotating the data sets used for calibration and prediction, until all four sets had been used for 
prediction. The results of the cross-validation are shown in Figure 9a. The prediction errors 
5 biases shown in Figure 9a are clustered near Omg/dl. Such clustering suggests that in the 
absence of a lamp change, this apparatus is capable of making satisfactory measurements of 
glucose concentration with these samples. 

Another cross-validation analysis was then performed. In this cross-validation 
analysis, the ninety-eight solutions discussed above were grouped into four subsets, and a 
n io different lamp was assigned for use as the illumination source for each subset. In this 
SJ analysis, data from three of the lamps was used to build a chemometric model to predict 

I : 1 

*jf glucose in data from the fourth lamp. This chemometric modeling procedure was repeated 

^ until each of the four data sets was used for prediction. The prediction results for the four 

q data sets are presented in Figure 9b. A comparison between the four data sets shows a very 

15 large lamp-to-lamp prediction bias. These results are again consistent with the findings 

u 

O presented in Figures 6 (the replacement of individual lamps) and Figure 7 (the modest 

rotation and/or translation of a single lamp by +/- 2 degrees), thus further illustrating the 
deleterious effects of interferents, such as illumination system variations, on the 
development of accurate chemometric models for preferred systems of the present 
20 invention. 

The apparatus diagrammed in Figure 10 is consistent with those discussed in 
detail with respect to Figure 5, with the clear identification of a radiation homogenizer 
90. In a preferred embodiment, the homogenizer 90 is positioned between the filter 78 



-30- 



hj 

ni 
\n 

; ; 5 



and the sample 40, as depicted in Figure 10. At this location, entering nearly 
monochromatic radiation is spatially and angularly homogenized prior to its distribution 
upon the sample 40. The placement of the homogenizer 90 at the above-described 
location is not to be construed as restricting the scope of the invention. The system 
5 depicted in Figure 10 is significantly simplified for illustrative purposes. Only certain 
specific elements within a far more elaborate spectroscopic system are diagrammed. All 
the elements depicted in Figure 10, however, are common to preferred spectroscopic 
systems of the present invention. The elements diagrammed, therefore, are to aid in 
identification of various aspects of the overall spectroscopic system. Thus, it should be 

10 understood that the present invention encompasses embodiments wherein various 
components of a spectroscopic system may be assembled in a relative order other than the 
one explicitly diagrammed in Figures 10 and 15. However, the homogenizer 90 is placed 
at a point between the emitter 14 and the tissue or sample 40, although other elements 
may be included between the homogenizer 90 and emitter 14 or between the 

15 homogenizer 90 and tissue or sample 40. This can also include the spectrometer 44, 
which in certain embodiments can be positioned between the emitter 14 and tissue 40. 

In a preferred embodiment, the radiation homogenizer 90 is a light pipe. Figures 
11a and lib show a perspective end view and a detail plan view of a light pipe 91 of the 
present invention. Light pipe 91 is generally fabricated from a metallic, glass 

20 (amorphous), crystalline, polymeric, or other similar material, or any combination 
thereof. Physically, the light pipe comprises a proximal end 92, a distal end 94, and a 
length 96 therebetween. The length of a light pipe 91, for this application, is measured 
by drawing a straight line from the proximal end 92 to the distal end 94 of the light pipe. 
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Thus, the same segment of light pipe 91 may have varying lengths depending upon the 
shape the segment forms. The length of the segment readily varies with the light pipe's 
intended application. 

In a preferred embodiment as illustrated in Figures 11a and lib, the segment 
forms an S-shaped light pipe. The S-shaped bend in the light pipe provides angular 
homogenization of the light as it passes through the light pipe. This conclusion is 
documented by the experiment and discussion associated with Figures 13a-c and 14a-c 
below. It is, however, recognized that angular homogenization can be achieved in other 
ways. A plurality of bends or a non-S-shaped bend could be used. Further, a straight 
light pipe could be used provided the interior surface of the light pipe included a diffusely 
reflective coating over at least a portion of the length. The coating provides angular 
homogenization as the light travels through the pipe. Alternatively, the interior surface of 
the light pipe can be modified to include dimples or "microstructures" such as micro- 
optical diffusers or lenses to accomplish angular homogenization. Finally, a ground glass 
diffuser could be used to provide some angular homogenization. 

The cross-section of the light pipe 91 may also form various shapes. In particular, 
the cross-section of the light pipe 91 is preferably polygonal in shape to provide spatial 
homogenization. Polygonal cross-sections include all polygonal forms having three to 
many sides. Certain polygonal cross-sections are proven to improve spatial 
homogenization of channeled radiation. For example, a light pipe possessing a hexagonal 
cross-section the entire length thereof provided improved spatial homogenization when 
compared to a light pipe with a cylindrical cross-section of the same length. 
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Additionally, cross-sections throughout the length of the light pipe may vary. As 
such, the shape and diameter of any cross-section at one point along the length of the 
light pipe may vary with a second cross-section taken at a second point along the same 
segment of pipe. 

In certain embodiments, the light pipe is of a hollow construction between the two 
ends. In these embodiments, at least one lumen or conduit may run the length of the light 
pipe. The lumens of hollow light pipes generally possess a reflective characteristic. This 
reflective characteristic aids in channeling radiation through the length of the light pipe 
so that the radiation may be emitted at the pipe's distal end. The inner diameter of the 
lumen may further possess either a smooth, diffuse or a textured surface. The surface 
characteristics of the reflective lumen or conduit aid in spatially and angularly 
homogenizing radiation as it passes through the length of the light pipe. 

In additional embodiments, the light pipe is of solid construction. The solid core 
could be cover plated, coated, or clad. Again, a solid construction light pipe generally 
provides for internal reflection. This internal reflection allows radiation entering the 
proximal end of the solid light pipe to be channeled through the length of the pipe. The 
channeled radiation may then be emitted out of the distal end of the pipe without 
significant loss of radiation intensity. An illustration of internal reflection and the 
resulting channeling is shown in Figure 12. 

Figure 12 depicts a plan view of a ray trace showing radiation 24 from a light 
source 14 (40-watt tungsten-halogen bulb) focused by an elliptical reflector 12 into, and 
through, a light pipe 91 of the present invention. In particular, Figure 12 illustrates how 
emitted radiation from a radiation source lamp is focused upon the proximal end of the 
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light pipe of the present invention. The focused radiation is internally reflected 
throughout the length of the light pipe. As the radiation is reflected, specific structural 
characteristics of the light pipe (here an S-shaped segment of hexagonal cross-sectioned 
pipe) angularly and spatially homogenizes the resulting radiation emitted at the pipe's 
distal end. 

Figures 13(a-c) are plots of the incidance of emitted radiation from the elliptical 
reflector and light pipe depicted in Figure 12. These plots have again been generated 
using TracePro V2.1, a ray trace program simulating the spatial distribution of emitted 
radiation from the radiation source emitter. More specifically, the plots of incidance are 
representative of the spatial distribution of emitted radiation at the distal end of the light 
pipe diagrammed in Figure 12. 

Figure 13a shows a plot of incidance of emitted radiation from the radiation 
source lamp coupled to the light pipe of the present invention. The resulting incidance 
plot is characterized by a substantial degree of spatial homogeneity. Spatial distribution 
of emitted radiation throughout the incidance plot varies slightly. A comparison of 
Figure 13a with that of Figure 3a illustrates the substantial improvement in spatial 
distribution throughout the incidance plot when using a light pipe of the present 
invention. 

Figure 13b shows a plot of incidance of the same radiation source lamp coupled to 
the light pipe of the present invention as depicted in Figure 13a, but after a 90-degree 
rotation of the filament producing the incidance plot. Again, the resulting incidance plot 
is characterized by a substantial degree of spatial homogeneity. In fact, there exist few 
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detectable differences in spatial distribution after the resulting 90-degree rotation as with 
the spatial distribution prior to the rotation. 

Figure 13c further depicts the spatial homogeneous distribution of emitted 
radiation using a light pipe of the present invention. Again, the spatial distribution in 
Figure 13c, after a one-millimeter translation, is very similar to those spatial distributions 
in Figures 13(a-b). 

Similar to Figures 13(a-c), Figures 14(a-c) show plots of the intensity of emitted 
radiation from the elliptical reflector and light pipe depicted in Figure 12. These plots 
have also been generated using TracePro V2.1 to simulate the angular distribution of 
emitted radiation from a radiation source emitter known in the art. More specifically, the 
plots of intensity are representative of the angular distribution of emitted radiation at the 
distal end of the light pipe diagrammed in Figure 12. 

Figure 14a shows a plot of intensity of emitted radiation from the radiation source 
lamp coupled to the light pipe of the present invention. The resulting intensity plot from 
the standard radiation source is characterized by a substantial degree of angular 
homogeneity. Angular distributions throughout the plot vary slightly. A comparison of 
Figure 14a with that of Figure 4a illustrates the substantial improvement in angular 
distribution throughout the intensity plot when using a light pipe of the present invention. 
For example, the "hole" in the center of the intensity plot caused by the glass nipple on 
the end of the radiation source lamp is no longer present and is now replaced with 
homogenized angular radiation. 

Rotating the filament of Figure 14a by 90-degrees produces an intensity plot 
illustrated by Figure 14b. Again, there are minor differences between the intensity plots 
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after, and prior to, the rotation. Translation of the filament of Figure 14a by one 
millimeter, as depicted in Figure 14c, once again documents reduction in variation of 
angular distribution as compared to the plots of Figures 14a-b. 

The ray trace plots of Figures 13(a-c) and 14(a-c) illustrate that the spatial and 
angular distribution of light at the output of the light pipe is highly stable with respect to 
modest translations and/or rotations of its filament. This is especially clear when comparing 
the ray trace plots of Figures 13(a-c) and Figures 14(a-c) using a light pipe of the present 
invention with Figures 3(a-c) and Figures 4(a-c) without the light pipe of the present 
invention. The light tube of the present invention has been effectively shown through these 
incidance and intensity plots to eliminate or substantially reduce the light source or 
illumination system as an interferent associated with chemometric modeling. It has been 
found that the use of the light pipe of the present invention allows construction of 
chemometric models of sufficient sensitivity to measure analyte concentrations. 

Another embodiment of the present invention is depicted schematically in Figure 15. 
In this embodiment, the tungsten halogen source 14 is placed at one focus of an elliptical 
reflector 1 10, and the proximal end 111 of a light pipe 1 12 is placed at the other focus 1 14. 
To improve the collection efficiency of the system, a separate back reflector 116 is 
positioned opposite the elliptical reflector 110 to capture and redirect light which would 
otherwise be lost from the system. The distal end 118 of the light pipe 112 then provides the 
source of radiation for the spectroscopic sample. 

Figures 16 and 17 show the simulated spatial and angular distributions of the light at 
the distal end 1 18 of the light pipe 112 of Figure 15. These distributions show substantially 
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improved homogenization as compared to the output of the standard system depicted in 
Figure 1. 

Another embodiment of the present invention is shown in Figure 18. In this 
embodiment, the tungsten halogen source 114 is placed at the focus 120 of a section of a 
parabolic reflector 122, and the proximal end 124 of a light pipe 126 is placed at the focus 
128 of a section of another parabolic reflector 130. The homogenized light exits the distal 
end 132 of the light pipe 126. The simulated spatial and angular distributions of the light at 
the distal end of the light pipe, shown in Figures 19 and 20, show substantially improved 
homogenization as compared to the output of the standard system depicted in Figure 1. 

Another embodiment of the present invention is shown in Figure 21. This 
embodiment is similar to the standard system depicted in Figure 1, except that the standard 
elliptical reflector has been replaced with a faceted reflector 140. This faceted reflector 140 
has the same general form as the elliptical reflector of Figure 1, but the smoothly varying 
shape of the standard elliptical form has been replaced with flat mirror facets 142 which 
locally approximate the standard shape. Such faceted reflectors 142 provide a high degree 
of spatial uniformity. Figure 22 is a simulated spatial distribution of the light at the second 
focus of the ellipse, showing substantially improved spatial homogeneity as compared to the 
output of the standard system of Figure 1. Figure 23 is a simulated angular distribution at 
the second focus of the ellipse which, unlike the other embodiments disclosed herein, 
exhibits a high degree of non-uniformity. 

The faceted elliptical reflector is an example of an embodiment of the present 
invention which produces only part of the desired characteristics in the output radiation. In 
the case of the faceted reflector 140, spatial homogenization is achieved but not angular 
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homogenization. In other cases, such as passing the output of the standard system through 
ground glass, angular homogenization is achieved but not spatial homogenization. In 
embodiments such as these, where only angular or spatial homogenization is produced (but 
not both) some improvement in the performance of the spectroscopic system may be 
expected. However, the degree of improvement would not be expected to be as great as for 
systems where spatial and angular homogenization of the radiation are simultaneously 
achieved. 

Another method for creating both angular and spatial homogenization is to use an 
integrating sphere in the illumination system. Although common to use an integrating 
sphere for detection of light, especially from samples that scatter light, integrating spheres 
have not been used as part of the illumination system when seeking to measure analytes 
noninvasively. In practice, radiation output from the emitter could be coupled into the 
integrating sphere with subsequent illumination of the tissue through an exit port. The 
emitter could also be located in the integrating sphere. An integrating sphere will result in 
exceptional angular and spatial homogenization but the efficiency of this system is 
significantly less than other embodiments previously specified. 

In order to evaluate the efficacy of the light tube of the present invention for 
reducing prediction error related to lamp variations, an experiment was conducted 
comparing a chemometric model using a light pipe of the present invention with a 
chemometric model without the light pipe of the present invention. The system of Figure 8 
depicts the system without the light pipe. Figure 24 is a diagrammed view of the system of 
the present invention for measuring glucose in scattering media having a tissue phantom 43 
as the sample source. The apparatus diagrammed in Figure 24 is consistent with that 
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discussed in detail with respect to Figure 8 except for the S-bend light pipe 91 which is 
included at the focus of the second silicon lens 74. 

The results of comparative testing between the system of Figure 8 and that of Figure 
24 which incorporates the light pipe are included in the box and whisker plots of Figures 25a 
5 through 25d. Figures 25a and 25b are duplicates of Figures 9a and 9b to provide easy 
comparison with the results included in Figures 25c and 25d. Thus, Figure 25a depicts the 
ability of the standard system with no bulb changes to predict glucose concentrations. 
Figure 25b depicts the system ability across four bulb changes. Figure 25c depicts the 
results of the system of Figure 24 across four bulb changes. Figure 25d shows the results of 
p 10 tests done on the system of Figure 24, but with the addition of a ground glass diffuser 78 
CQ prior to the light pipe 91 . Figures 25c and 25d clearly show that the embodiments of Figure 

24 are highly effective in improving the predictive accuracy of the apparatus and 
chemometric model over the system of Figure 8. Further, the greatest benefit is derived 
when the ground glass diffuser 78 and the S-bend light pipe 91 are used together which 
15 results in the highest degree of homogenization of the light incident on the sample. 

The performance of the illumination system to a known radiation emitter difference 
can be quantified. A method for quantifying the performance of the illumination system is 
to create both angular and spatial distribution plots under two known but different 
conditions. The differences between the two similar metric plots can be quantified. The 
20 known emitter difference to be used for quantification is a one-millimeter translation of the 
lamp filament. 

Angular and spatial distribution plots can be created by using standard ray trace 
packages such as TracePro V2.1 or through direct measurements. The image of the 
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illumination system beam can be measured by using any standard intensity mapping scheme 
and by using a goniometer. This allows both the spatial and angular distributions of the 
illumination output to be determined. 

Optical modeling or direct measurement of the system should occur before and after 
movement of the filament. In order to standardize the calculation for many applications, the 
image should be divided into approximately one hundred equally sized "bins" (or squares), 
with ten bins across the diameter of the output image. This requirement is easily satisfied 
when performing ray trace analysis and can be accomplished by either measuring the beam 
in a ten by ten grid or by sampling at finer spacing and then averaging the data. The spatial 
and angular distributions for the initial emitter state are then subtracted from the 
corresponding distributions after movement of the lamp filament by one millimeter. The 
resulting images represent either the angular or spatial variance that occurred due to the 
emitter perturbation. In order to quantify the angular or spatial variance, all the data in the 
different images are put into a vector for easier calculation, and the vector is normalized so 
that its length equals 1 . This normalization is achieved by dividing each data point by the 2- 
norm fll-lb), which is equivalent to the Euclidean distance of the vector, 



where X is the vector of the difference image and n is the number of data points in 
the vector. 

The normalization step ensures that the magnitude of every difference-image is 
comparable. Following the normalization step, the standard deviation of the normalized 
image vector is calculated, and this metric is an indication of the amount of variance 
introduced by the known emitter difference, 
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The standard deviation of the normalized image vector for both angular and spatial 
distributions was calculated for three different illumination systems. 

1. Acceptable System: This illumination system is a light source (40-watt 
tungsten-halogen bulb) focused by an elliptical reflector into a ground 
glass diffiiser, specified as a weak angular homogenizer, with subsequent 
coupling into a hexagonal light pipe with a length to diameter aspect ratio 
of 3 to 1. The system is modeled such that the filament image fully fills 
the input into the hexagonal light pipe. 

2. Preferred System: the illumination system is the same as the acceptable 
except that the length to diameter aspect ratio is 7 to 1. 

3. Ideal System: The illumination system is composed of a light source (40- 
watt tungsten-halogen bulb) focused by an elliptical reflector into a ground 
glass diffuser, specified as a strong angular homogenizer, with subsequent 
coupling into an s-bend hexagonal light pipe with a length to diameter 
aspect ratio of 7 to L The system is modeled such that the filament image 
fully fills the input into the hexagonal light pipe. 

Based upon testing with these three illumination systems, the degree of 
homogenization can be generally classified as acceptable, preferred and ideal. Table 2 
shows the standard deviations of the spatial distribution for the three systems. Table 3 
shows the standard deviation for angular distribution. 
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TABLE 2 



Vertical Filament Filament Rotation 
Acceptable 0.053 0.050 

Preferred 0.045 0.042 

Ideal 0.039 0.034 



TABLE 3 



Acceptable 

Preferred 

Ideal 



Vertical Filament 
0.044 

0.032 

0.027 



Filament Rotation 
0.066 

0.054 

0.050 



There is another metric that can be used to evaluate the efficacy of an illumination 
system in reducing error inflation following bulb changes. This metric is known as the 
multivariate signal to noise (mSNR). The typical signal to noise (S/N) calculation is a 
univariate measure; it is defined as the maximum signal in a spectrum divided by the 
standard deviation of the baseline noise. 

When a multivariate calibration is used, the signal from two or more wavelengths 
is used to quantify the analyte of interest. Because of this, unless the noise is random or 
'white' noise, the standard deviation of the baseline (as used in univariate S/N 
calculations) is an inexact inappropriate noise estimate. Furthermore, the maximum 
signal in the spectrum is also an inexact inappropriate measure of the overall signal since 
the multivariate calibration uses signals from multiple wavelengths. The mSNR metric, 
however, uses the multivariate net analyte signal and the error covariance matrix and 
therefore gives a better estimate of the signal to noise for multivariate calibrations. 
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The net analyte signal is that part of the analyte spectrum which is orthogonal 
(contravariant) to the spectra of all interferents in the sample. If there are no interfering 
species, the net analyte spectrum is equal to the analyte spectrum. If interfering species with 
similar spectra to the analyte are present, the net analyte signal will be small. Because the 
calibration depends on the net analyte signal, the multivariate signal to noise metric takes 
this measure into account. 

The mSNR can be calculated if two pieces of information are known. The net 
analyte signal (NAS) for the analyte of interest must be known, but this may be estimated 
from the regression vector, b (the model), 

NAS = A- (3) 

ihi; 

where ||.||2 represents the 2-norm of the vector. 
The error covariance matrix (E), which describes the error structure of the multi- 
wavelength spectral data, is also needed for the mSNR calculation, 

E = s T *s (4) 
where 6 is a vector containing the noise at each wavelength. 

x = x 0 + e (5) 
where x is a measured spectrum, x 0 is the "true" spectrum in the absence of noise, 
and e is the noise. 

The error covariance matrix, E, measures how noise is correlated across 
wavelengths. The spectra used to calculate the error covariance matrix are spectra that 
have a constant amount of the analyte of interest and are obtained or processed in a 
manner to identify the spectral variances due to the variance of interest. In practice, a 
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repeat sample should be used and the only variance introduced into the system should be 
the spectral variance being identified. In this invention, the variance source of interest is 
spectral variances due to emitter changes. Thus, spectral data from a repeat sample is 
obtained under different emitters. If the noise is uncorrelated, the error covariance matrix 
will have no off-diagonal elements, but in many cases, this will not be true. In such 
cases, the error may 'overlap' spectrally with the net analyte signal. In other words, this 
will introduce 'noise' into the measurement of this particular analyte. The 'Noise' may 
be calculated as, 



Noise = Vv T Sv (6) 



10 where 



n NAS 

yj v = 



|NAS| 2 

The mSNR at unit concentration may then be calculated by, 



(7) 



llNASlI llNASlI 

mSNR = * ^= ", " 2 (8) 

Noise Vv T ]£v 

The inverse of the net analyte spectrum, llmSNR, is an estimate of how much error will 
15 be added to prediction estimates if the type of noise in e is present in the spectra being 

used to predict the analyte concentration (or other property). 

When an illumination system is insensitive to emitter variances, there will be little 

effect on the spectral noise; in other words, the error covariance matrix, E, will be close 

to diagonal. In that case, the mSNR will be high. In the case where the system is 
20 sensitive to emitter variances or source fluctuations, correlated noise will be introduced 

and that will create off-diagonal elements which will be present in the error covariance 
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matrix E. When these spectral variances or noise interferes (co-varyies) with the net 
analyte signal, the mSNR gets smaller and its inverse increases. 

Table 4 shows the mSNR and 1/mSNR values calculated for four different 
illumination systems. These systems include a standard system with no bulb changes, the 
preferred embodiment system (with s-bend light pipe and diffuser) and also one that 
contained a straight light pipe (acceptable system). 

TABLE 4 

System mSNR 1/mSNR 

No bulb change (Ideal level) 0.2 5 

Bent light pipe & diffuser (Preferred level) 0.033 30 

Straight light pipe only (Acceptable level) 0.0166 60 



CD It is clear that bulb changes influence each system differently. The mSNR is 

10 highest when no bulb change occurs, and lowest when the standard system with limited 
source homogenization is used. Conversely, the greatest inflation in prediction errors can 
be seen in that system (approximated by 1/mSNR). 

These mSNR values were calculated using the study measuring the 98-solution set 
that was described previously. The NAS was calculated using the model (b) generated 
15 from the data set where a single bulb was used (equation (1)). This model had no 
knowledge of bulb changes, and so the net analyte signal corresponds to that in the 
absence of source fluctuations. For each illumination system, there were four bulb 
changes as described before. For each bulb, in addition to the 90-solution set, additional 
'repeat' samples were measured. These 'repeats' were simply samples that contained all 
20 of the analytes at concentrations at the center of the calibration. Thus, to isolate the 
spectral variance due to bulb changes the spectral data was processed in the following 
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manner. Multiple 'repeat' spectra at each bulb were measured, and the average repeat 
spectrum for each bulb was calculated using these data, hereafter referred to as the 
average bulb spectrum. Each average bulb spectrum can be thought of as the 'x' in 
equation 5. The mean repeat spectrum is simply the average spectrum of the average 
bulb spectra. To calculate the error, 6, associated with each bulb, the mean repeat 
spectrum was subtracted from the average bulb repeat spectra, 

n 

I*. 

£,=X,-J=!— (9) 

n 

where n is the number of bulbs in the analysis (4 in this example). The S matrix was then 
calculated using equation 4, and equations 6-8 were then calculated to find the mSNR. 

Now referring to Figure 26, another aspect of the present invention is depicted. The 
system depicted provides spectral filtering or bandpass filtering to eliminate unnecessary 
wavelengths or bands of wavelengths from the light prior to contact with the tissue. This is 
accomplished by placing one or more elements between the light source and tissue. The 
elements can include absorptive filters fabricated of any transparent or partially transparent 
substrate; single layer or multi-layer dielectric coatings deposited on any transparent or 
partially transparent substrate; a grating or prism which disperses the radiation, permitting 
unwanted wavelengths to be blocked from reaching the tissue; and/or an aperture which 
selectively blocks undesirable radiation. 

A preferred system for bandpass filtering is depicted in Figure 26 which depicts a 
light source 100 placed within an electrical reflector 102. Figure 26 also depicts a 
hexagonal S-bend light pipe 104 to receive light from the source 100. A series of filters are 
placed between the light source 100 and the light pipe 104. The first optical filter is a silicon 
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filter 106 which is anti-reflection coated to transmit at least ninety percent (90%) of the in 
band incident light. The silicon filter passes wavelengths of light longer than 1.1 micron. 
The second optical filter is preferably a KOPP 4-67 colored glass filter 108 that, in 
combination with the silicon filter, passes light in the 1.2 to 2.5 micron spectral region. The 
slope of the KOPP filter is such that is preferentially passes light at wavelengths between 2.0 
and 2.5 micron. The third optical filter is an ORIEL WG295 absorption filter 1 10 that cuts 
out wavelengths longer than 2.5 micron. The front surface of the WG295 filter can be 
polished or finely ground. If the front surface is finely ground, the WG295 acts as a dififuser 
as well as a light filter. It has been found that these filters prevent burning of the tissue, 
while enhancing the signal-to-noise ratio of the system by band limiting the light to only the 
spectral region of interest. The effect of band limiting the light is to reduce shot noise 
generated by the photon flux incident on the detector. 

An alternative combination of filters to achieve spectral bandpass filtering is 
depicted in Figure 5. With this embodiment, the two silicon lenses 72,74 absorb 
wavelengths shorter than approximately 1.2 microns and longer than approximately 10 
microns. The cyan filter 76 is an absorptive filter such as a Hoya CM-500 to absorb mid- 
infrared radiation at wavelengths of approximately 2.8 microns and longer. Further, a 
SCHOTT WG-295 absorptive filter 78 is included to absorb radiation at wavelengths 
approximately 2.7 micron and higher. Figure 27 graphically depicts the individual and 
combined spectral transmission of the components shown in Figure 5, along with the 
"spectral fingerprint" of glucose. As depicted in the graphs, this combination of absorptive 
filters and silicon lenses acts to block unwanted wavelengths, while still permitting 
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transmission of radiation in the glucose fingerprint region. Similar combinations of filters 
can be utilized based on analytes of interest to be analyzed. 

It is also recognized that other modifications can be mad^efthe present disclosed 
system to accomplish desired homogenization of light. Fpf'example, the light source could 
be placed inside the light pipe in a sealed arrangement which would eliminate the need for 
the reflector. Further, the light pipe^ould^ replaced by an integrator, wherein the source is 
placed within the integrator as disclosed in U.S. Patent Application Serial No. , 



entitled "Encoded Multiplex Variable Filter Spectrometer," filed on the same date herewith 
and incorporated^ reference. Further, the present system could be used in non-infrared 
10 application^ achieve similar results in different wavelength regions depending upon the 
type (^analysis to be conducted. 

Having thus described the preferred embodiments of the present invention, those of 
skill in the art will readily appreciate that yet other embodiments may be made and used 
within the scope of the claims hereto attached. 
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